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Abstract 

The linearized Vlasov equation for a plasma system in an external 
constant magnetic field and the corresponding linear Vlasov opera- 
tor are studied. The solution of the Vlasov equation is found by the 
resolvent method. The spectrum and eigenfunctions of the Vlasov op- 
erator are also found. The spectrum of this operator consists of two 
parts: one is continuous and real; the other is discrete and complex. 
Interestingly, the real eigenvalues are uncountably infinitely degener- 
ate, which causes difficulty solving this initial value problem by using 
the conventional eigenfunction expansion method. It also breaks the 
natural relation between the eigenfunctions and the resolvent solution 
that the eigenfunctions can be considered as the coefficients of e~' lujt 
in the Laplace (or resolvent) solution. 
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1 Introduction 



In plasma physics, many interesting phenomena such as plasma oscillations, 
instabilities and Landau dampings can be studied and understood through 
a very simple model of plasma. In this model, ions are assumed to be mo- 
tionless and form a positively charged background and the collisions between 
electrons are neglected. In 1945, Vlasov proposed a nonlinear equation, now 
known as Vlasov equation, to describe this model |]]. Since then, a lot of 
effort has been made to study this equation. When there is no external 
field, one usual way to treat this problem is to first linearize the equation 
by assuming the system is very close to the equilibrium state, then reduce it 
to a one-dimensional equation. Many interesting phenomena can be under- 
stood through this simple one- dimensional linearized Vlasov equation, fj, |3|, |J 

In this paper, I shall study this system in an external constant magnetic 
field. After the Vlasov equation is linearized, two methods, the eigenfunction 
expansion method and the resolvent method, are tried to solve the equation. 
The resolvent method is proved to be successful. The eigenfunction expan- 
sion method is, surprisingly, not successful even after all the eigenf unctions 
of the corresponding linear operator, Vlasov operator, are worked out. The 
difficulty is caused by the fact that the real eigenvalues are uncountably in- 
finitely degenerate. 

It has been studied by Case f| and Arthur et al || the explicit relation 
between the Laplace transform (or the resolvent) approach and the eigen- 
function expansion approach to the one-dimensional Vlasov equation with 
no external field. Especially, the latter demonstrated how to construct the 
van Kampen-Case modes || [7| through the resolvent solution. Their results 
show that the van Kampen-Case modes are, loosely speaking, just the coef- 
ficients of e luJt in the resolvent (or Laplace) solution. This is no longer true 
when we don't reduce the Vlasov equation to a one-dimensional problem or 
simply we can't because of the existence of external fields. The eigenmodes 
(or eigenfunctions) corresponding to the real eigenvalues presented in Section 
3 can not be constructed with the resolvent solution in Section 5, which is 
caused by the infinite degeneracy mentioned above. 

The plan of this paper is as follows. In Section 2, the derivation of 
the linearized Vlasov equations and its related linear Vlasov operators are 
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sketched, primarily to introduce notations. The eigenfunctions are presented 
and discussed in Section 3 without detailed calculations. Next, the related 
adjoint problem is discussed briefly and the orthogonal relation between the 
original eigenfunctions and the adjoint eigenfunctions is proved. In Section 
5, the exact solution is given by using the resolvent method. 



2 Vlasov Equations and Vlasov Operators 

This simple model of plasma is described by the Vlasov equation 
d d F d 

— p(x, V, t) = -V- T^p(x, V,t) -^ZP(X, V,t) (1) 

ot ox m e Ov 

where p(x, v, t) is the single distribution function of electrons and F is the 
force felt by a single electron. In this paper, the system is in an external 
constant magnetic field B , so 

F = -ev x B -eE (2) 

-» Of p(x' V* t) 

where E(x, t) = e— / dx'di? - - — ^-j— is the electric field at position x and 
ox J \x — x'\ 

time t generated by the electron system itself. If the system is just slightly 
away from an equilibrium state fo(v) 

p(x, v, t) = n f (v) + f(x, v, t) (3) 

where no is the average density and / nofo, then it can be described very 
well by the linearized Vlasov equation which is 

-^ = - v -J^ + —( vxB °yj^ + — / dx dv ■ -J* ( 4 ) 



dt dx m e dv m e dx J \x — x' 

Taking the Fourier transform 

f(x, v ,t) = fk(v,t) exp (ifc • x ) (5) 

we have for the A;-component function f%(v, t) 

[ ]f t fk(v, t)= k- vf u (v, t) + i^-(v x B ) ■ ^f u (v, t) 



(6) 



where tu p = (A7re 2 n /m e ) 1 ^ 2 is the plasma frequency. To simplify the notation, 

— * 

let us drop the index k and let 

then equation (6) becomes 

i-f^,t) = k-vf(v,t) + i^(vxB,)-—f(v,t)- V (v) Jtiff{tf,t) (8) 

This is the equation I shall study mostly in this paper. Consequently, the 
corresponding linear Vlasov operator /C is defined as 

K,f{v) = k ■ vf(v) + i^-(v x B ) ■ —f(v) - V (v) J dv'fiv') (9) 

This is an integro-differential, unbounded and non-self-adjoint linear opera- 
tor. In terms of /C, equation (8) can be put into a compact form 

i^mt)=icmt) (io) 

The equation for the case where B = is 

i^mt)=Kof(v,t) (11) 

with 

1Cof{v) = k ■ vf(v) - V (v) J dv'fiv') (12) 

— * 

Integrating over the two components of v perpendicular to k, we can reduce 
(11) to a one-dimensional problem 

i^/M)=£o/M) (13) 

With _ _ roc 

fCofiv) = kvf(v) - fj(v) / dv'f(v') (14) 

J— oo 

— * 

Here v is the component of v along the direction of k and k is the magnitude 
of k. The bar over functions indicates that they have been integrated over 
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the two components of v perpendicular to k. This convention is assumed 
throughout this paper. Equation (13) is the one which has been studied 
thoroughly by van Kampen, Case and many others 0, [7], §. 

Formally, equations (10), (11) and (13) are just the same as the Schrodinger 
equation. But the Hamiltonian operator appearing in the Schrodinger equa- 
tion is self-adjoint and bounded from the below while fco, /Co and /C are 
unbounded and not self-adjoint. These are all typical equations arising in 
linear evolution systems. There are many standard methods to study this 
class of linear equations. For example, equation (13) can be solved by the 
Laplace transform |2], || , the resolvent method and the eigenfunction ex- 
pansion method [[5], [TJ. I shall try to use the eigenfunction expansion method 
and the resolvent method to solve equation (10) and (11). 

3 Eigenfunctions and Eigenvalues 

In this section, I shall present the eigenfunctions and spectra of operators 
JCo, /Co and /C. These three operators are similar in many respects as they 
are supposed to be. Their spectra are the same: continuous real eigenval- 
ues, discrete real eigenvalues and discrete complex eigenvalues. Since the 
existence of discrete real eigenvalues depends on the choice of the equilib- 
rium function fo(v), I shall only consider, for simplicity, the case where the 
discrete real eigenvalues don't exist. Also, their eigenfunctions are similar: 
most of them are singular containing 8 functions. However, there is one big 
difference: the real eigenvalues of /Co and /C are infinitely degenerate while 
the eigenvalues of JC are not degenerate. This infinite degeneracy causes the 
difficulty expanding functions in terms of these eigenfunctions and relating 
these eigenfunctions with the resolvent solutions. 

Operator KLq: 

The eigenfunctions of K,q were first constructed by van Kampen |F] then com- 
pleted by Case || and are now widely known as van Kampen-Case modes. 
Their results can be summarized as follows. 
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The eigenequation for eigenvalue z is 

/oo 
dv'g z (v') = zg z (v) (15) 
-oo 

It is linear so the functions can be normalized as 

/oo 
dvg z (v) = 1 (16) 
-oo 

When the eigenvalue is real, say, it is v>, the eigenfunction is 

g v [y) = V-^~ + \(v)5(kv - V) (17) 
KV — V 

Here V means the principal value integration. The normalization condition 
(16) requires 

The discrete complex eigenvalues Vj (j = 1,2, ...,m ) are determined by 

e (z) = 1 - r dv-^± = 1 _ [ m J® = o (19) 

7-oo kv-Z J k-V-Z 

The corresponding eigenfunction is 

which satisfies normalization condition (16). Eigenf unctions (17) and (20) 
are the famous van Kampen-Case modes. It is clear that all the eigenvalues 
of /Co are not degenerate when the van Kampen-Case modes are normalized 
according to (16). 

Operator KLq: 

The eigenequation is 

fco9z(v) = k-vg z {v) - 7](v) j di?g z (i?) = zg z (v) (21) 

As this operator is linear, we set the normalization condition 

J dvg z (v) = 1 (22) 
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The eigenfunction corresponding to a real eigenvalue v is 

g „(0) = V - ^ + X(v)5(k -v-v) (23) 
k ■ v — v 

The normalization condition (22) requires that X(v) satisfy 

/ dv\(v)5(k ■ v - v) = 1 - / dvvj^— (24) 
J J k ■ v — v 



It is very clear here that there are infinitely many choices for X(v) to satisfy 
the above condition. For example, all functions X(v) = a(vi,v 2 )\(v) satisfy 
(24) as long as 

/'OO roo 
dv 1 / dv 2 a(v 1 ,v 2 ) = 1 (25) 
-oo J —oo 

where V\ and v 2 are the two components of v perpendicular to k. Obviously, 
the choices of a(vi,v 2 ) are uncountably infinite, which means there are in- 
finitely many corresponding eigenfunctions for any real eigenvalue v. In other 
words, the degeneracy of every real eigenvalue is uncountably infinite. 

The complex eigenvalues of /Co are the same as JCo, Vj [j = 1,2, ...,m ), 
the zeroes of e (z), whose corresponding eigenfunctions are 

9iW=9» i W = T 3 ®— (26) 
k ■ v — Vj 



Operator /C: 

I shall present the results directly without going into the detailed calculations. 
The method of computing eigenfunctions is very similar to the one given in 
the Appendix. First, let us set up a coordinate system and introduce some 
notations: 

B = B z , k = k±x + k\\z (27) 
v = v± cos9x + v± sin 9y + v \\z (28) 
Thus the operator 1C becomes 

Kf{v) = (k ± v ± cos9 + k m )f(v) - wo^m ~ V(v) J difftf) (29) 
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where c^o = eB /m e is the cyclotron frequency. The same normalization 
condition 

J dvG z (v) = 1 (30) 

is set for the eigenequation 

K,G z {v) = zG z (v) (31) 

As the complex eigenvalues of /C are determined by 6q(z) = 0, the complex 
eigenvalues Zj (j = 1,2, ...,m) of /C are determined by e(z) = where 

e(z) = 1 + vr £ r v ± dv ± f dr,,, + ^ W + 2 ^H (32 ) 

KJ ^ h J-oo 11 2(* - mwo - Vll) 

Here and after, all summations are assumed to be over all integers. J n 's are 
the Bessel functions J n (k±v±/cuo). Two functions, rj± and r)\\, are defined as 

V± = %k ± £-f {if) , m = ^±f (if) (33) 

Note the equilibrium state function fo(v) has been assumed to be a function 
of only v± and t> y , and m, the number of the discrete complex eigenvalues, is 
possible to be infinite. The eigenfunction corresponding to Zj is 

G 3 {v) = eMi k±V±Sin6 ) E e -™ {Jn+1 + Jn - l)r]± { 2 + J ^ (34) 
For a real eigenvalue /i, we have 



G,(v) = exp(i^^) £ e'^ Jn+1 + ^ W2 + J ^ 



U n At — nCJ - fc||D|| 

/- \ .,k ± v ± Bm6 
2 y an(W||, f±) exp — 1( nu)d{j2 — nooo — k\\ v \\) 

n UJ 



(35) 



where arbitrary functions as long as they satisfy the only constraint, 

the normalization condition (30). This means that there are infinitely many 
eigenfunctions G^{y) corresponding to each real eigenvalue \i. The infinite 
degeneracy arises again. It is easy to check that all these expressions reduce 

— * 

to the ones for B Q = . 
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We see that there are many similarities among K.q, /Co and /C. But there 
is one important difference: the real eigenvalues of K,q are not degenerate 
while the real eigenvalues of both /Co and /C are uncountably infinitely de- 
generate. As it has been shown by Case ||, the solution of equation (13) can 
be expanded in terms of the van Kampen-Case modes. In other words, we 
can use the eigenfunction expansion method to solve equation (13) just as 
we often do to solve the Schrodinger equation. However, this method is not 
suitable for solving equations (10) and (11) even though we know the eigen- 
functions of /Co and /C. Due to the infinite degeneracy, there is no obvious 
way that functions can be expanded in terms of these eigenf unctions. It is a 
common belief that for an initial value problem of a linear evolution system, 
to find the solution is equivalent to find the eigenfunctions and spectrum of 
the linear operator in this system. Seemingly, this is not the case for /Co and 
K. 



4 Adjoint Equations and Operators 

Operator K is not self-adjoint, so it is interesting to know the functions 
orthogonal to its eigenfunctions. For this purpose, we consider the adjoint 
equation of (10) 

i^ffrt) = KfW (36) 
where the linear operator /C is defined as 

Kj{v) = k ■ vf(v) + i—(vx B ) ■ —f{v) - [ d^(t7W) (37) 

The discrete complex eigenvalues of K are also determined by equation (32). 
Its eigenfunction corresponding to Zj is 

r< t-f\ c k±v± sin 6, ^ _ in6 J n . . 

G>) = eXP(l ~^ } ? e ^--r^o-Vll { ) 

while its eigenfunction of corresponding to a real eigenvalue fi is 



u V nuJ o - h v \\ 

(39) 

+ 2_^a n (v\\,v±) exp — 1( n&)d{fi — nu — k\\V\\) 
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where subject only to the normalization condition 

J dv V (v)G^v) = 1 (40) 

It means that the real eigenvalues of K. are also infinitely degenerate. There- 
fore, we see that the adjoint operator K has the same spectrum structure 
and similar eigenf unctions as the original operator JC. By straightforward 
substitution, it can be proved 

J dv[jCf(v)]*g(v) = J dvf*(v)[ICg(v)} (41) 

where "*" represents the complex conjugate. This relation leads easily to the 
orthogonal relation between the eigenf unctions of /C and K 

{z* - z') J dvG* z (v)G z ,{v) = (42) 

One interesting point to note is that for the complex eigenvalues we have 

JdvG* Z] (v)G Zj (v)=0 (43) 

since z* ^ Zj. The adjoint equation of (11) can be considered as a special 
case of B = 0, no additional treatment is necessary. 



5 Resolvent Method 

As it has been pointed out in the previous sections, the conventional eigen- 
function expansion method is not suitable for solving equation (10). We have 
to resort to other methods. It turns out that the resolvent method is a suc- 
cessful choice as it is for the one dimensional equation (13) 0. The resolvent 
of /C is defined as 

n*) = -\ (44) 

where z is a complex variable. If the initial function is f(v, 0) then the 
function at time t is 

m *) = ^ / dze- izt K(z)f(v, 0) (45) 

Z7T1 J 
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where the integration contour surrounds all the singularities of 1Z(z). For 
most operators, how their resolvents act on a function can be found only 
approximately. In this case, the explicit form of resolvent is found exactly ( 
the detailed derivation is given in the Appendix ), which is 

*M/(9, 0) = F(,, 0) - G(,, (46) 



where 



k±v±sm6 ^ _ ine (J n +i + J n -i)ri±/2 + J n 
G(z,v) = exp(i )z^ e 



G(z,v) = exp(i )J2 



uj „ -2 — ncj — k\\V\\ 



.k ± v ± sm0 x ^ _ in6 J n 



u ^ z-nuj - k\\v\\ 



F(z, v) = exp(i^^) £ f m J n e~^° 



with 



z-(m + n)cu - fc||U|| 



Similarly, for /C we have 



1 ^ Jdggj)^ 
-~f(v,0) = F(z,v) -G(z,v) (47) 



z-JC JK ' ' v ' 7 v ' ' e(z) 

As is well-known, the singularities of the resolvent TZ(z) give the spectrum 
of /C. It is not hard to see from (46) that the singularities of TZ(z) include 
the zeroes of e(z), discrete real poles, nuo + k\\v\\, and a branch cut along the 
real axis. Therefore the spectrum of K is just what we obtained earlier: all 
real values and a set of discrete complex values determined by e(z) = 0. Let 
us write down the solution explicitly by choosing a specific contour 

(48) 

+ 2^iL dZe i F M-G(z,v) — } 
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where C± run parallel to the real axis and are chosen such that all the singu- 
larities of K.(z) are enclosed between them. We can use the residue theorem 
to calculate the above integration. If we are only interested in the solution 
for t > 0, then two big semi-circles can be attached from below to C± to make 
two closed contours 0±. Since there are no singularities inside the contour 
0_ the second part of the integration (48) is identical to zero. In contrast, 
+ encloses singularities: poles, such as Zj, and a branch cut along the real 
axis. To get rid of the troubling branch cut, we can replace the functions 
in the integrand with their plus analytic continuations ||. These plus con- 
tinuations are identical to their original functions in the upper-half complex 
plane but are different in the lower-half plane. They have no branch cuts 
but have new poles different from the original poles in the lower-half plane. 
These new poles are the mathematical origin of the Landau dampings 0. [| 
in this plasma model. 

When B = 0, equation (48) becomes 

1 f(v 0) n(v) I dv^B 

/(iT, 0) = J{ 1 } - V[ J f=p- (49) 



z — fco z — k-v z — k-v e o( z ) 

Therefore, the solution of equation (11) is 

mt)= -±[ ize -^ { mx- m /d ' 7gs } 



2iriJc + ' z-k-v z-k-v e Q (z) 

+ ± f dze -i*t { MO) _ V(v) Sdv^ 
z-k-v z-k-v e (z) 



(50) 



After the integration over the two components of v perpendicular to k, the 
above equation can be reduce to one dimensional 



/<„,*)= -±[ ^-"{im- m /d "^ } 



2tti Jc+ z — kv z — kv €o(z) 

2wi Jc- z — kv z — kv €o(z) 

which is the solution of equation (13). 



(51) 
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It is important to note here that the eigenfunctions g z (v) of ICq can be 
constructed through its resolvent solution (51). The way was demonstrated 
clearly in ||: Contour C± is divided into two parts, small circles surrounding 
poles and two straight lines approaching to the real axis from both below and 
above. The first part gives easily the eigenfunctions gi{v) corresponding to 
the complex eigenvalues. The second part leads to the real eigenfunctions 
g v {v). For operators /C and )C, their complex eigenfunctions gj(v) and Gj(v) 
can be constructed easily similarly through their resolvent solutions (48) 
and (50), respectively. But straightforward calculation immediately shows 
the technical difficulties for the real eigenfunctions g^v) and G^{v) be con- 
structed out of (48) and (50) in a similar way. It is not worthwhile to dwell 
on the technical difficulties which involve very complicated calculations and 
formulas since one may argue there may be another method to circumvent 
these difficulties. Let us focus on the underlying reason, the uncountably 
infinite degeneracy. Usually, the eigenfunctions of a linear operator can be 
labelled exactly either by a set of real numbers or integers or both. For op- 
erators /Co and JC, I still labelled their eigenfunctions by real numbers and 
integers. But the labelling is not exact since there are infinite number of 
eigenfunctions for each real number. It certainly requires the introduction of 
non-trivial measure and other concepts in functional analysis to expand any 
function in the vector space expanded by the eigenfunctions g^{v) or G^{v). 
All of this is unlikely to be achieved by some clever algebraic manipulations 
with (48) and (50). 



6 Conclusion 

In summary, I have studied the linearized Vlasov equation of plasma in an 
external constant magnetic field Bq and found its exact solution by the re- 
solvent method. Also the corresponding linear Vlasov operator K. has been 
discussed and its eigenfunctions are found. The spectrum of this unbounded 
operator has two parts: continuous real eigenvalues and discrete complex 
eigenvalues. This is very similar to the spectrum of Kq, the one dimensional 
linear operator derived from the linearized Vlasov equation without external 
fields |7|, |!|. But there is a striking new feature for K, that the real eigen- 
values of K, are uncountably infinitely degenerate. This new feature leads to 
the difficulty in expanding functions in terms of the eigenfunctions of /C. 
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Appendix 

In this appendix, the detailed derivation leading to (46) is given. The eigen- 
functions Gj(v) and G^(v) can be found in the similar way. The method 
is called the integration along the characteristic curves || or the integra- 
tion along the "unperturbed orbits" ||]. In the following calculations, it is 
assumed that lm.z > 0. To make formula compact, let us set two notations 

F={z-K)f{v) (Al) 

E = exp(ifc • x - izt) (A2) 
It is easy to verify with the aid of (9) that 

4 + » ■ I - x *> ■ = fe - E "I f < A3) 

where r\ = r](v), f = f(v) and the non-specified integration is over velocity. 
Along the "unperturbed orbit", i.e., the orbit in the phase space of one 
electron moving in magnetic field Bq, we have 

dt dt " dx rrig" °^ dv ^ 
Then equation (A. 3) becomes 

i^(fE) = FE-r ] E J f (A5) 

Integrating it, we get 

if = f dt'F'E'E- 1 - ( f dt'ri'E'E- 1 ) f f (A6) 

J—oo J— oo J 
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where the primed functions mean the functions of primed variables. Note in 
the above integration over time t, the fact Imz > has been used to choose 
— oo as the lower limit. Another integration with respect to velocity gives 

f f _ II^dt'F'E'E-^ 

J i + ff^dt'rj'E'E-i V ■ } 

Plugging this back into (A. 6), we have 

{z-K)- l F = -i f dt'F'E'E' 1 

J —oo 

, i IJ'dt'F'E'E' 1 



(A8) 



+ / dt'rj'E'E- 1 - 
J-oo i 



Adopting the coordinate system and notations in (27) and (28), we write out 
(A.8) explicitly 

f-OO 

K(z)f(v) = -tj o drf{v h v ± ,e-oo t)e-^ 

+ r drnme-^ 1 d€S ~ dr/( "" l ^ ° - ^ 
Jo aTJ l^> e l-ijdvj^dr^e-^ 



where 

0(r) = -^( s in0 - sin(# - u t)) + (k\\v\\ - z)r (A 10) 

= v± c °s(^ - uot) + v\\ ( A - n ) 

The case Imz < can be dealt with similarly, which also leads to (A. 9). The 
step from (A. 9) to (46) is direct integration with respect to r using Bessel 
functions. 
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